# Required functions: `rsyrs.mean`

set.seed(0)
N=400; n=2; k=200; id=1:N
periodic=sin(id/3.18)+.01*rnorm(N)
s1=seq(sample(1:k,1),N,k); s2=seq(sample(1:k,1),N,k)
s3=seq(sample(1:k,1),N,k); s4=seq(sample(1:k,1),N,k)
s5=seq(sample(1:k,1),N,k); s6=seq(sample(1:k,1),N,k)
s7=seq(sample(1:k,1),N,k); s8=seq(sample(1:k,1),N,k)
s9=seq(sample(1:k,1),N,k); s10=seq(sample(1:k,1),N,k)
p1=periodic[s1]; p2=periodic[s2]; p3=periodic[s3]
p4=periodic[s4]; p5=periodic[s5]; p6=periodic[s6]
p7=periodic[s7]; p8=periodic[s8]; p9=periodic[s9]
p10=periodic[s10]
rsys=data.frame(p1,p2,p3,p4,p5,p6,p7,p8,p9,p10)
rsys.means=apply(rsys,2,mean)

muhat.rsys('periodic',rsys.means,k,N)